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Abstract We present an inflationary model preceded by a bounce in a metric f(R) 
theory. In this model, modified gravity affects only the early stages of the universe. 
We analyse the predicted spectrum of the gravitational waves in this scenario using 
the method of the Bogoliubov coefficients. We show that there are distinctive (os- 
cillatory) signals on the spectrum for very low frequencies; i.e., corresponding to 
modes that are currently entering the horizon. 



1 Introduction 

We propose a bouncing scenario 0] within the context of a metric f(R) theory 0: 
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S=— 2 jd\^-gf{R)+SW. (1) 

The bounce we will consider in our model is followed by an inflationary era which is 
asymptotically de Sitter where, in addition, the gravitational action approaches the 
Hilbert-Einstein action on that regime, such that the modification to Einstein's Gen- 
eral Relativity (GR) affects exclusively the very early universe, around the bounce 
and a few e-folds after that. We will constrain the model obtaining the spectrum of 
the stochastic gravitational fossil as would be measured today. 
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2 Model for the early universe 



Inspired on the de Sitter solution for a closed space-time, we define the scale factor 
around the bounce as: 

a(t) = flfe cosh (Hi n ft ) , (2) 

where a(t) is the scale factor, a/, is a constant quantifying the size of the universe at 
the bounce. The parameter H m f is related to the energy scale of inflation just after 
the bounce. 

In a Friedmann-Lemaitre-Robertson- Walker (FLRW) universe with a spatially 
flat metric and the scale factor defined as in Eq. (|2}, the minimization of the f(R) 
action (fTJ leads to a second order differential equation for the function /. Solving 
this equation in conjunction with appropriate physical constraints gives |3|: 



/(r)=2// i 2 nf v / r^3cos 



\fi> ( 9-r\ /3r 
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In the above equation r=Rj H? nf . 



n — arccos — arccos \ / 

V 3 / V2r-3 



(3) 




RIHm RWm 

Fig. 1 These plots show: (left) the behaviour of f(R) as a function of R/H? nf (see the blue curve) 
and linear approximation R — 6H? nf , well during inflation (see the red curve); (right) the behaviour 
of fg (see the continuous curve) and fgg (see the dashed curve) as functions of R/ H? nf . 



3 Energy spectrum of the gravitational waves 

The spectrum of the gravitational waves is determined using the method of the con- 
tinuous Bogoliubov coefficients a and /3, as in Ref. H. The graviton density of the 
universe is given by |j3| 2 , while the dimensionless logarithmic energy spectrum of 
the gravitational waves (GW) of angular frequency co is defined at the present time 
T7o as Q: 



The spectrum of gravitational waves in an f{R) model with a bounce. 



3 



o 




s 
o 
G 



-16 -15 

Log 10 w (rad/s) 




-16 -15 -14 
Log lo ai (rad/s) 






-10 






-12 






-14 




o 
G o 


-16 






-18 




O 
J 


-20 






-22 






-24 




3 




17 




Log 10 w (rad/s) 



-15 
Log 10 m (rad/s) 



Fig. 2 These plots show the spectra of the energy density of GW at present time: (top left) Ob- 
tained with an f(R) treatment for the perturbations (blue continuous curve) and a GR treatment for 
the perturbations (red dashed curve). Ejnf = 1.5 x 10 16 GeV. fj n j = — arccosh(10)// i ^ nf 1 . a b = 2 x 10 3 . 
(top right) Otained in an f(R) treatment for different values of a b . The value of a b increases 
from the red curve to the blue curve as: a b = 2x 10 2 ; a b = 2 x 10 3 ; ci/, = 2x 10 4 ; a b = 10 x 10 5 . 
a b = 10 x 10 6 . Em = 1.5 x 10 16 GeV. / ini = -arccosh(10)H i ; f 1 . (bottom left) Obtained in an f(R) 
treatment for different values of t m \. The value of t\„\ increases in absolute value from the red curve 



to the blue curve as: ?j, 



-arccosh(10)// inf ' ; 1\ r 
il6r 



-ai-ccosh(20)// inf l ; t\ 



-arccosh(50)//j 



inf • 



f ini = -arccosh(100)//r/ . E mt = 1.5 x 10 K, GeV. a b = 2 x 10 3 . (bottom right) Obtained in an/(R) 
treatment for different values of £; n f. A comparison is made between the results obtained with a 
fixed value of a\, (continuous curves) and a fixed value of ai,H m \ (discontinuous curves). The value 
of £jnf increases from the red curve to the blue curve: E m f = 1.5 x 10 14 GeV; £; n f = 0.5 x 10 15 GeV; 
Em = 1.5 x 10 15 GeV; E inf = 0.5 x 10 16 GeV; E inf = 1.5 x 10 16 GeV. (aw = 10a 6 . a b = 2 x 10 3 ). 



To calculate the evolution of the gravitational waves, we express the the contin- 
uous Bogoliubov coefficients in terms of the variables X and Y, see Refs. 0121, 
which obey the set of differential equations: 

X" = (k 2 - — ^ X, and X' = -ikX. (5) 

Here, k is the wave-number, a prime indicates a derivative with respect to the confor- 
mal time 7] (a = dt /drj) and z = <1\[Jr- The differential equations (0 are integrated 
from an initial time f; n ;, set before the bounce, until the present time. We describe 
the late time evolution of the universe in a GR setup, using the ACDM model 
complemented with a radiation phase and making the connection between the f(R) 
driven early inflation and the radiation phase with a model of a modified General- 
ized Chaplygin Gas suitable for the early universe J7|. The results obtained for the 
GW spectra are shown in Fig. [2] 
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The existence of the bounce in the early universe affects the spectrum of GWs only 
in the low frequency range, < l(T 12 Hz, where various peaks appear whose position 
and intensity depend on the parameters of the mode. The fact that the oscillatory 
structure appears in the spectra of (i) the GR treatment and (ii) the f(R) treatment 
suggests it is not a consequence of the effects of / (R) -gravity. Similar oscillations 
have been obtained in works of loop quantum cosmology first pointed out by Afonso 
et al (9J- Due to the low energy density of the cosmological GWs, the results ob- 
tained in this work are hard to be detected in the near future (see Fig. 2 of Ref. IflOl 
and Fig. 6 of Ref. ifTTI ). The detection of the B-mode polarization of the CMB ra- 
diation seems to be the best candidate to obtain information about the cosmological 
GWs. 02) 
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